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We consider a one-dimensional, time-reversal-invariant system with attractive interactions and 
spin-orbit coupling. Such a system is gapless due to the strong quantum fluctuations of the su¬ 
perconducting order parameter. However, we show that a sharply defined topological phase with 
protected, exponentially localized edge states exists. If one of the spin components is conserved, 
the protection of the edge modes can be understood as a consequence of the presence of a spin 
gap. In the more general case, the localization of the edge states arises from a gap to single particle 
excitations in the bulk. We consider specific microscopic models and demonstrate both analytically 
and numerically (using density matrix renormalization group calculations) that they can support 
the topologically non-trivial phase. 


I. INTRODUCTION 

Topological phases of matter are typically character¬ 
ized by a gapped bulk spectrum, and protected gapless 
edge states with unique properties. The existence of a 
finite energy gap in the bulk plays a crucial role in the 
topological protection of the edge states. This therefore 
raises the question whether a topological phase, i.e. a 
phase with protected exponentially localized edges states, 
can exist in a gapless system [1]. 

Superconductors have been shown to host a variety 
of topological phases [2-7], depending on the symme¬ 
tries and the dimensionality of the system. In one spa¬ 
tial dimension, such a phase hosts protected edge modes, 
termed Majorana bound states, which are in particular 
interesting due to their non-Abelian exchange statistics. 
Realization of this phase requires proximity coupling 
a one-dimensional system to a bulk three-dimensional 
superconductor [8-12]. In truly one-dimensional sys¬ 
tems with intrinsic attractive interactions, strong quan¬ 
tum fluctuations of the superconducting order parameter 
leave the system gapless. E.g., this situation can arise if 
a quantum wire is coupled to a superconducting wire. 
It was shown that topological protection is then much 
weaker in general [13, 14] (the edge states are generically 
only power-law localized, rather than exponentially local¬ 
ized) and depends on the microscopies of the system [15]. 

In the presence of time-reversal symmetry a different 
kind of topological superconductors can be realized [16- 
26]. If the bulk is fully gapped, these phases host a 
Kramers pair of Majorana bound states [24, 27-29]. 

In this work, we show that one-dimensional time rever¬ 
sal symmetric systems can support a well-defined topo¬ 
logical superconducting phase, even when the particle 
number is conserved. The two necessary ingredients to 
realize this phase are attractive effective interactions and 
spin-orbit coupling. We start by giving a field theoreti¬ 
cal argument for the existence of the topological phase, 
distinct from a trivial one. We show that the boundary 
between the topological and the trivial phases hosts ex¬ 
ponentially localized edge states and discuss their prop¬ 


erties. The edges exhibit an anomalous relation between 
time reversal and the local fermion parity operator, just 
like in the fully gapped case. We then demonstrate the 
existence of the topological phase and its properties nu¬ 
merically. To this end we consider a simple model that 
exhibits a phase transition between the trivial and the 
topological phases and study it using the density matrix 
renormalization group (DMRG) technique. Finally, we 
consider a system consisting of a semiconducting wire 
with spin-orbit coupling and repulsive interactions cou¬ 
pled to a superconducting wire, and show using bosoniza- 
tion and weak coupling renormalization group (RG) that 
it can be driven into the gapless topological phase. 

II. FIELD THEORY OF GAPLESS 
TIME-REVERSAL-INVARIANT TOPOLOGICAL 
SUPERCONDUCTORS 

A. Spin conserving case 

We begin by considering an interacting, time rever¬ 
sal invariant one-dimensional electron gas (IDEG), de¬ 
scribed at low energies as a Luttinger liquid. In bosonized 
language, the Hamiltonian can be written as 

^0 = E ^ (1) 

where p, a correspond to the charge and spin degrees of 
freedom, p{x) = —^dx4>p and s^(x) = —i^dx4>a are the 
charge and spin densities, respectively, and ^9a is the 
field conjugate to (fa, [4>a {x) ,da {x')] = iirQ^x' — x), 
where 0(x) is the Heaviside step function. 

The relation to fermionic operators, describing modes 
linearized around the Fermi momentum kp, is given by 

R{L)^ _ ,^-i{±^<j>p-ep+s{±^4>,,-e„)) ^ 2 ) 

where s denotes the spin of the fermion. Us are the Klein 
factors that impose the anti-commutation relations be¬ 
tween the different spin species, -|- (—) signs correspond 
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to R(L), representing right(left) movers, respectively, and 
a is the short distance cutoff in the theory [30]. Back- 
scattering processes give rise to cosine terms that can 
gap out some of the modes in the system. We assume 
the system to be at a generic filling, such that there are 
no relevant umklapp processes. 

Our system is symmetric under time-reversal (TR) 
symmetry, denoted by T (class Dill in the Zirnbauer- 
Altland classification [31], such that = —1). Un¬ 
der time reversal, —)■ L^, —)■ —R^, —)■ i?^, 

and Ri —)■ —L^. These relations correspond to the fol¬ 
lowing transformations of the bosonic fields: 9p —>■ —9p, 
'Pp ^o- 9cr, (per ~4’cn U'f and Ui —> —U'f. 

These transformation rules reproduce the correct behav¬ 
ior of the spin and charge densities and currents under 
TR. 

We consider a system with spin-orbit coupling, and 
hence no SU(2) spin symmetry. Let us assume, for sim¬ 
plicity, that the spin is conserved along one direction (e.g. 
the z direction, such that is conserved). This condi¬ 
tion will be relaxed later on. The most general four- 
fermion back-scattering interaction consistent with time- 
reversal symmetry is 

g (^R\l^lIr^ + h.c)j = cos (2(/)^), (3) 

where the coupling g is a real number. Note that such a 
cosine term is invariant under T, according to the trans¬ 
formation rule of (per above. (Higher order processes are 
also possible, but are less relevant in the RG sense.) If 
this cosine term is relevant, it opens a spin gap, driving 
the system into a “Luther-Emery phase” [32] with only 
one gapless (charge) mode. The sign of g determines the 
nature of the resulting phase. To understand this re¬ 
call that the spin-singlet and spin-triplet pairing order 
parameters are given by 

Oss = R\lI + l\rI (X cos {cp ^), 

O^S = R\lI - l\rI oc sin {cP ^), 

while the spin-density wave and charge-density wave or¬ 
der parameters are given by 

OcDW = oc e“*‘^'’cos {(pc) , 

= r\L^ - oc e-*‘^'>sin {(Pc ). 

For g < 0, the cosine term pins the field (pc to zero (or 
equivalently to any integer multiple of tt, i.e. tto with 
n S Z), resulting in the dominant superconducting cor¬ 
relations being of the spin-singlet order parameter. (As 
the charge sector remains gapless, these pairing correla¬ 
tions decay with a power law dictated by the Luttinger 
parameter Kp and no true long range order can develop). 
For g > 0, (pc is pinned to tt (n + i), where n S Z, and 
the dominant superconducting correlations are now of the 
spin-triplet order parameter. These two cases correspond 


to two distinct phases; in order to go between them with¬ 
out breaking time-reversal symmetry, the coupling g has 
to cross zero, resulting in the closing of the spin gap [33]. 
We refer to these phases as trivial and topological respec¬ 
tively, in analogy with the fully gapped case. When the 
superconducting order parameter is conventional s-wave, 
we expect the system to be in the trivial phase, while for 
a p-wave order parameter the system is topological. This 
identification is in agreement with the vacuum being in 
the trivial phase, as a large back-scattering potential at 
the end of the system pins (pc to zero. 

The phase diagram of an interacting IDEG with spin- 
orbit coupling, and the possibility of a phase with dom¬ 
inant triplet superconducting correlations was discussed 
in Refs. [34, 35], where it was demonstrated how this 
phase can arise in a model that describes certain quasi- 
one-dimensional organic conductors. However, to the 
best of our knowledge, the topological nature of this 
phase (manifested in its protected edge modes, as we 
argue below) has not been discussed. 

Consider an edge of a topological system, or equiv¬ 
alently a boundary between a trivial and a topological 
phase. Since the field (pc is pinned to 7rni on one side 
of the boundary and to 7r(n2 -I- ^) on the other, where 
ni ^2 are integers, there must be a kink of minimal mag¬ 
nitude (pc across the boundary. Such a kink in (pc 

corresponds to an accumulation of spin: 

■Sz = y s^{x)dx = -J ^dx(pcdx = ±^, (6) 

i.e. half the spin of an electron. The same “fractional 
spin” appears at the edge of a time-reversal invariant 
fully gapped topological superconductor [24]. 

The lowest energy configurations of (pc for a topological 
system of finite size are shown in Fig. 1. The two config¬ 
urations plotted in blue solid lines correspond to an even 
number of particles in the system with total (S'z) = 0. 
These configurations are degenerate up to a splitting ex¬ 
ponential in system size, as in the case of a fully gapped 
topological superconductor with time-reversal symmetry. 
Adding a single particle in the bulk of the system requires 
creating a kink of magnitude ±7r (as the extra particles 
carries spin Sz = i^)- Such a kink costs a finite amount 
of energy. Adding a particle near the edge, however, only 
requires flipping the direction of the ^ kink at that edge 
without an extra energy cost. Since the number of par¬ 
ticles in the system is a good quantum number, adding 
a particle costs charging energy, but this contribution 
decreases with the system size L as The resulting 
configurations with an odd number of particles and to¬ 
tal (Sz) = are given by red dashed lines in Fig. 1. 
(Note that the two states with an odd number of par¬ 
ticles are exactly degenerate, in accord with Kramers’ 
theorem.) We therefore expect that the gapless topolog¬ 
ical phase will host single particle edge excitations which 
cost zero energy in the thermodynamic limit. This is in 
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FIG. 1: Low energy configurations of (f>cr in a finite size gap¬ 
less topological superconductor. Configurations plotted in red 
dashed lines correspond to total spin of ±| in the system, i.e. 
odd number of particles, while configurations plotted in blue 
solid lines correspond to zero total spin, i.e. even number of 
particles. To calculate the local tunneling density of states 
we calculate the matrix element of a single particle creation 
operator (x) between the states denoted by |0) and |1). 



(a) (b) 


FIG. 2: The ground state energy of a system with open 
boundary conditions as the number of particles is varied for 
(a) a system in the trivial phase (b) a system in the topo¬ 
logical phase. In both cases the contribution of the chemical 
potential is subtracted, A indicates the single particle gap. In 
the topological phase a low energy state appears for an odd 
number of particles due to the presence of the low energy edge 
state. Note that a single particle gap Is still present in the 
bulk of the system. The data is obtained using DMRG study 
of the model given by Eq. 16 for a system of length L = 100 
sites with model parameters t = 1, U = —1 and spin-orbit 
coupling w = 1 in (a) and w = 0 in (b). The chemical potential 
that is subtracted in both cases is found from a linear fit of 
Ejv to N. The red solid lines are fits to parabolas. The finite 
curvature of the parabolas is due to the charging energy, and 
depends on the size of the system. 


sharp contrast to the gapless trivial phase where no such 
edge excitations exist. 

This distinction can be seen in Fig. 2, depicting the 
ground state energy of a system with open boundary con¬ 
ditions as a function of the number of particles in the two 
cases. For a system in the trivial phase, states with an 
odd number of particles lie on a parabola separated by A, 
the single particle gap in the system, from the states with 
an even number of particles. For a system in the topo¬ 
logical phase, an extra particle can be added at the edge 
of the system at a cost of the charging energy only. The 
two states are distinguished by the pair binding energy 
in the system 


Eb = E 21 V +1 — 2 (-l^2N+2 + E2 n) (7) 


as the size of the system is varied. Extrapolating to the 
infinite system size limit, we expect to obtain the single 
particle gap, limi_>oo7?B = A > 0, if the system is in the 
trivial phase. If the system is in the topological phase and 
with open boundary conditions, we expect WniB^aoEB to 
be zero. In contrast, in a system with periodic boundary 
conditions, both phases are characterized by a finite Eb 
in the thermodynamic limit. 

Experimentally, the existence of a low energy single 
particle edge state is refiected in the local tunneling den¬ 


sity of states (TDOS) 

—ImG^ = 

TT 

(|(0|^'J (x) |n)f-b Kola's (a;) |n)|^) , 

n,s 

( 8 ) 


in the w —>■ 0 limit. Here |0) denotes the ground state 
of the system and |n) an excited state at energy Hun- 
Assuming the number of particles in the ground state is 
even, the configuration of in this state is given by one 
of the blue solid curves in Fig. 1, e.g. the one denoted 
by |0). The lowest energy excited state that contributes 
to the sum in Eq. 8 corresponds to a configuration of 
(jiu denoted by 11). Using the bosonized representation of 
dfl {x) and a mode expansion for the bosonic fields, (pp^n 
and 0p,o-, we obtain for the matrix element: 


(lld/j ix)\0) 



a N Try 


\ 

l) 

Vx/ 


( 9 ) 


where a = + \Kp, is the spin gap, and ^ ^ 

is the correlation length in the spin sector. For details of 
the calculation see Appendix A. 


B. Non-spin conserving case 

We now generalize the analysis of the previous section 
to the more generic case, in which none of the spin com¬ 
ponents are conserved. 
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FIG. 3: Generic spectrum for a system without inversion sym¬ 
metry. i?(L)i ,2 denote the right(left) movers in each band. 
Under time-reversal symmetry Ri —> L 2 and R 2 —>■ —Li. 

Let us first consider a case where the system has no 
inversion center. Then, the single particle spectrum is 
generically of the form shown in Fig. 3. We linearize 
the modes close to Fermi energy and denote by R{L)i ,2 
the right (left) movers in each band. Under time-reversal 
symmetry Ri —)■ L 2 and i ?2 —t —Li. We can therefore 
assign a pseudo-spin f, | to the two bands, respectively, 
and write the linearized fermionic modes in terms of the 
bosonic fields (j)p^cr and exactly as before, where 4>cr 
is now related to the density of the pseudo-spin. The 
only four-particle back-scattering term allowed by time- 
reversal symmetry is R\l^l\R 2 + h.c., or cos(2(()o-) in 
terms of the bosonic fields, and all the analysis presented 
for the Sz conserving case still holds. 

If the system has an inversion center, the resulting 
phase diagram is richer. The case of spin anisotropic 
interactions in an inversion symmetric system with time- 
reversal symmetry was analyzed in [35] . Here, we outline 
this analysis briefly, and comment on the topological na¬ 
ture of the various phases and their boundary properties. 

After an appropriate rotation in spin space, the most 
general four-fermion inversion symmetric interaction can 
be written as 

iLint = A, + ^ (10) 

d—x,y,z 

where A^ = (fs^) R is a singlet order parame¬ 

ter, Aj = L'^ (^is^d • s^ i? is a time-reversal symmetric 
triplet order parameter characterized by the orientation 
of a headless d vector, and we have dropped terms that 
do not couple the right and left movers. Here L{R) are 
spinors in spin space, L{R)'^ = {L{R)^, L{R)i). Note 
that inversion symmetry prohibits cross terms that cou¬ 
ple the singlet and triplet order parameters. 

The back-scattering part of this interaction can be 
written in terms of the bosonic fields as 

- 9s) cos {2(l)a) + {gy - gx) cos {49^}} . 

( 11 ) 

The resulting phase diagram hosts four phases with dom¬ 
inant superconducting correlations. These correspond to 
pinning of either (j)^ or da with two possible values for 
each, depending on the sign of [g^ — gs) or (gy — g^), re¬ 
spectively. Let us denote the phase in which = 0 
by SS (dominant singlet superconducting correlations), 
and the phase with (j)„ = hy TSz (dominant triplet 


correlations with d || z). The other two phases with 
9c = n'^,{n + 4) will be denoted as TSx and TSy, 
respectively. 

The SS and TS^ phases are exactly the ones discussed 
in the context of an Sz conserving system. These phases 
are distinct from each other (in the sense that they can¬ 
not be adiabatically connected without a phase transi¬ 
tion in which the gap in the spin sector closes) as long 
as time reversal symmetry is preserved. Below, we argue 
that the TSx,y,z can be adiabatically connected to each 
other without breaking time reversal symmetry. How¬ 
ever, they are distinct as long as certain mirror symme¬ 
tries, Mx^y^z (where is a mirror symmetry that takes 
X —> —a;, etc.) are maintained; in the presence of such 
a mirror symmetry, an interface between a pair of TS 
phases with different spin quantization axes hosts pro¬ 
tected edge states. 

To see that the TSx and TSy can be adiabatically 
connected without TR symmetry breaking, consider the 
following family of interaction Hamiltonians: = 

— \g'\ AjAj, where d = cos{a)x + sin(a)y. In bosonized 

form, oc — |g'| cos (4(0^, — q;/ 2)); therefore, as 

a varies from 0 to the value to which 9c is pinned 
changes from 0 to J, the values that correspond to the 
TSx and TSy phases, respectively. Interpolating from the 
TSx or TSy phases to the TSz phase is more difficult in 
the context of Abelian bosonization; nevertheless, by a 
rotation in spin space we can show that these phases can 
also be adiabatically connected in a similar fashion. 

Note, however, that this path in Hamiltonian space de¬ 
fined above includes terms like At A-, that break AL, and 
My. In the presence of either of these mirror symmetries, 
such terms are forbidden, and the TSx and TSy phases 
are distinct from each other. 

In order to understand the topological nature of the 
TSx and TSy phases, we examine the low energy spec¬ 
trum in a system with open boundary conditions. Let us 
consider a system of length L, with boundary conditions 
imposed by a strong back-scattering potential at both 
ends. The back-scattering term has the form 

-\Vb\ 'y' RIL^+ h.c. = - cos (0p) cos {(j)c) ■ (12) 
^' 7ra 

Such a potential pins both (j)p and (j)c to either an even 
or an odd multiple of tt at both ends of the system, a: = 0 
and L. A system described by a state with 9c pinned to a 
specific value in the bulk can not satisfy these boundary 
conditions, due to the non-trivial commutation relations 
[(j)c (0), 9c (a;)] = in, where 0 < x < L. 

Moreover, the fermion parity V = , where N = 

— i {(j)p (L) — (j)p (0)) is the total number of particles, is 
not well defined in such a state. To see this note that 
using the boundary conditions dictated by Eq. 12, we 
can rewrite V in terms of the spin degrees of freedom as 
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since {(j)p (L) - (pp (0)) = {(p^ [L) - (pa (0)) mod 27r. 

Therefore, the following anti-commutation relations hold 

j}giea{x) _ 

Ground states of the, e.g. TSx phase, with well defined 
fermion parity can be constructed in the following way: 
1 ^ 0 ,±) = \0a = 0) ± lOa = Tt) and I'i/'f ,±) = \^a = f) ± 
\9a = ^), where the symmetric (anti-symmetric) su¬ 
perpositions correspond to even (odd) fermion parity 
states. (Recall that all the states [do- = n|) are de¬ 
generate in the bulk.) Note also that these two sets of 
states are distinct and are related by a transformation 
that tranfers a single fermion between the two edges, 
'I'l (L) (0) ^ e± 5 [ 0 <T(i)± 0 <T(o)] ^ This is reminiscent 

of the four-fold ground state degeneracy present in the 
fully gapped time-reversal-symmetric topological phase, 
linked to the existence of a low energy single particle ex¬ 
citation at the edge of the system. In our gapless system, 
however, states with a different total fermion parity also 
differ in their particle number, and hence are not degen¬ 
erate; a system with a given particle number has only 
two ground states, given e.g. by |V'o,+) and \%p^^+). 


where xi _2 are arbitrary points within the bulk of the 
system, a few correlation lengths away from the left and 
right edge, respectively. We have used the fact that 
^g-*( 0 <T(a; 2 )- 0 a( 2 :i))^ _ since (pa is pinned to a constant 
value in the bulk. The phase of the parity operators was 
chosen such that — 1- (Recall that at the 

boundary of the system (pa is pinned to an integer multi¬ 
ple of TT, while in the bulk it is pinned to tt (n -I- |) ; hence 
(pa (xi) - (pa (0) = (pa (L) - (pa {X2) = f mod TT.) Since 
the time-reversal operator T is anti-unitary and takes (pa 
to —(pa, we obtain {T,Vr^l} = 0 . 

For a system in the trivial phase, a consistent phase 
choice for the parity operators that ensures = 1 

is 

giving the usual commutation relations [T,= 0. 

III. NUMERICAL EVIDENCE FOR THE 
EXISTENCE OF A TOPOLOGICAL PHASE 


C. Symmetry fractionalization at the edges 


To better understand the nature of the topological 
phase, we consider the connection between time rever¬ 
sal and the local fermion parity at its edges. In fully 
gapped systems, symmetry-protected topological phases 
are characterized by a projective (fractionalized) rep¬ 
resentation of the symmetry operators when acting on 
the low-energy states of the edge [36-41]. In the case 
of a time-reversal symmetric topological superconductor, 
there is an anomalous relation between time reversal and 
fermion parity at the edge: TVr^l = —'Pr^lT, where 
Pr.l is the local fermion parity operator acting on the 
right or left ends of the system, respectively [16, 42]. 
These are local operators acting near the two edges, 
defined such that in the low-energy subspace the total 
fermion parity operator has the form V = Vr'Pr- We 
will show that, even though our system is gapless, the 
gap to single fermion excitations in the bulk guarantees 
that a similar decomposition of fermion parity in terms 
of local edge operators holds; therefore, the topological 
phase is characterized by the same anomalous relation 
between time reversal and fermion parity as in the mean- 
field case. 

For concreteness, we demonstrate this for a system in 
the TSz phase. In the low energy subspace, the total 
fermion parity given by Eq. 13 can be written as a prod¬ 
uct of local fermion parities at the edges of the system 
with 

(14) 

— _l^-^(.rj>„{L)-<p^{x2)) ^ ^ 


We next demonstrate the existence of the topological 
phase in a simple model, using the density matrix renor¬ 
malization group (DMRG) [43-46]. 

Gonsider spinful electrons on a ID lattice, with spin- 
orbit coupling in an alternating direction and spin- 
anisotropic interactions that explicitly favor triplet pair¬ 
ing 


H — Hq + Hsoc + Hint 


Ho 
Hsoc = iv 


h.c. 


I \ 


i-K-i)" 


i-(-iF 


+h.c. 


-) 


s,s' j + 


Hint — ^3 — 

( 16 ) 

where t is the hopping amplitude along the chain, v is 
the spin-orbit coupling strength and U is the interaction 
strength. 

In the absence of the single-particle spin-orbit cou¬ 
pling, i.e. for u = 0 , the z component of the spin is 
conserved. For any t/ < 0 a spin gap opens and we ex¬ 
pect to find the system in the topological phase. The 
spin-orbit coupling term is chosen such that it breaks 
spin symmetry completely. We will use it to study the 
robustness of the topological phase without spin conser¬ 
vation, and its range of stability. 

We calculated the low energy spectrum for a system of 
length L = 100 with N = 20 particles and open boundary 
conditions. The parameters used in this calculation are 
t = 1, U = —1, u = 0. The ground state is two-fold de¬ 
generate, as expected in the topological phase. We could 
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FIG. 4: DMRG results for the model described by Eq. 16 
in absence of spin-orbit coupling, u = 0. A system of length 
L = 100 sites with model parameters t=l, ?7 = —lis 
considered, (a) Expectation value of the z component of the 
spin along the chain in the ground states of the system for even 
and odd number of particles. The blue and green solid curves 
correspond to the two degenerate ground states of a system 
with = 20 particles and zero net spin. Accumulation of 
spin at the edge of the system is observed, with the integrated 
spin in the left (right) half of the system being ±|. The red 
dashed curve corresponds to the ground state of a system 
with A = 21 particles and total spin S' = -|-|. Due to the 
non-zero spin localized at the edge, and power-law decaying 
spin density wave correlations expected in a phase with (pa- 
pinned to ^ (see Eq. 5), a spin density wave pattern is formed 
in the bulk, (b) Matrix elements for the transition between 
each of the two ground states with A = 20 particles and 
the ground state with A = 21 particles by adding a spin up 
particle at a position i along the lattice. As can be seen, the 
matrix elements are non-zero only at either end of the system 
depending on the initial state. This is in agreement with the 
existence of low energy single particle states at the edges of 
the system. 


not resolve the expected exponential splitting between 
the two lowest energy states; this is presumably because 
the correlation length is much smaller than the system 
size. The two states found by DMRG are the minimally 
entangled states with an integrated spin of (S'z) = 
near the two edges. The configuration of the z com¬ 
ponent of the spin in the two ground states is shown in 
Fig. 4a. The topological phase is expected to have power- 
law decaying spin density wave correlations in the bulk 
(see Eq. 5 for pinned to |). As a result, the spin 
polarization at the edge induces a spin density wave that 
decays as a power law into the system, clearly visible in 
the figure. 


We then obtain the lowest energy state of the sys¬ 
tem with an extra spin up particle. Let us denote it 
by |2A + 1). We calculate the matrix elements between 
each of the states with an even number of particles, |2A), 
to the state |2A -I- 1) by adding a spin up particle at a 


position i along the lattice, 


(2A+l|cI_,^|2A) 


(see Fig. 


4b). We find that the matrix elements are non-zero only 
at either end of the system, depending on the initial state. 
This is in agreement with the existence of a low energy 



FIG. 5: Pair binding energy (see Eq. 7) as function of 
system size for the model Hamiltonian given in Eq. 16 with 
parameters t = 1, U = —I, spin-orbit conplings v = 0 and 
u = 1, and fixed density n = ^ = 0.2. The red solid lines 
are fits to parabolic curves. For v = 0 the system is in the 
topological phase and the pair binding energy tends to zero 
as 1/L —>■ 0. For v = 1 the system is in the trivial phase and 
Eb tends to a finite constant equal to the single particle gap 
in the system. 


single particle state at the edge of the system. 

In addition we calculated the pair binding energy (de¬ 
fined in Eq. 7 above), as function of system size, keeping 
the density of particles fixed at n = 0.2 (see Eig. 5). We 
find that it indeed tends to zero as the system becomes 
large, as expected in the topological phase. This once 
again indicates the existence of a low energy edge state, 
with an energy going to zero in the thermodynamic limit. 

For non-zero u, Sz is no longer conserved. More¬ 
over, conservation of Sz (or any other spin component) 
mod 2 is also broken, ensuring no residual symmetries 
are present. However, we find that the system remains in 
the topological phase for a finite range of spin-orbit cou¬ 
pling strength, v < Vc- To see this, we perform finite-size 
scaling of the energy gap (defined as the energy differ¬ 
ence between the first excited state and the ground state 
in an even particle number sector); keeping the density 
constant, we choose system sizes for which the number of 
particles is even. We find the gap to be exponentially de¬ 
creasing with system size, as expected in the topological 
phase. For each spin-orbit coupling strength, we extract 
the inverse correlation length in the system by fitting 
the energy gap, AE, vs. system size, L, to the form 
AE = We find that the correlation length ^ 

diverges as v approaches Vc « 0.22 (see Fig. 6). 

For V > Vc we calculate the pair binding energy, Eb, 
for different system sizes. As can be seen from Fig. 6, the 
binding energy tends to a non-zero value that increases 
with V, as the system size is increased (see also Fig. 5 for 
Eb as function of system size for u = 1). This indicates 
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FIG. 6: Phase diagram obtained for the model described by 
Eq. 16 with model parameters t = 1, U = —1 and total par¬ 
ticle density n = —0.2. As the spin-orbit coupling 

strength is increased, a phase transition from the topologi¬ 
cal to the trivial phase is observed. For v < Vc ~ 0.22 the 
gap (defined as the energy difference between the first excited 
state and the ground state) in a system with an even number 
of particles decreases exponentially with the system size, as 
expected for the topological phase. Data points marked by 
a blue cross correspond to the inverse correlation length ob¬ 
tained from the finite size scaling of the gap AE = 

For V > Vc we plot the the pair binding energy, Eb (see Eq. 
7), for different system sizes. In the limit L —>■ oo the pair 
binding energy tends to a non-zero value, increasing with v, 
indicating an opening of a trivial single particle gap in the 
system. 


the opening of a single particle gap and the absence of 
edge states in the trivial phase, as discussed in the pre¬ 
vious section. 

To emphasize the difference between the trivial and the 
topological phases, we calculate the ground state energy 
of the system, as the number of particles is varied in 
each phase (see Fig. 2). For v = 1 > Vc, the system 
is in the trivial phase with odd particle number states 
having energy larger by A, the single particle gap, with 
respect to the states with an even number of particles. 
For V = 0 < Vc the system is in the topological phase 
and there is no gap for single particle excitations; the 
extra particle can be accommodated at low energy at the 
edges. 


IV. POSSIBLE REALIZATION OF THE 
TOPOLOGICAL PHASE - QUANTUM WIRE 
COUPLED TO A ID SUPERCONDUCTOR 


Finally, we discuss a possible realization of the gap¬ 
less topological phase in a composite semiconducting¬ 
superconducting one-dimensional system. 

Consider a semiconducting quantum wire with re- 



FIG. 7: (a) Realization of the topological phase in a semicon¬ 
ducting quantum wire with spin-orbit coupling and repulsive 
short-range interactions, U > 0, coupled to a one-dimensional 
superconductor with intrinsic short-range attractive interac¬ 
tions V < 0. Spin-orbit coupling is assumed to have similar 
magnitude and opposite sign on the two wires, and the chem¬ 
ical potential is taken to be such that the inner modes in 
both wires are at fc = 0 to simplify the analysis, (b) Energy 
dispersions in the two wires and the relevant and marginal 
scattering processes. Two independent pair tunnelings of the 
inner and the outer modes are denoted by gi and g 2 respec¬ 
tively. The later is suppressed due to the mismatch in the 
energy dispersions for each of the spin flavors, as a result of 
the spin-orbit coupling. Back-scattering process due to the 
repulsive interactions in the semiconducting wire is denoted 

by g- 


pulsive electron-electron interactions coupled to a one¬ 
dimensional s-wave superconductor, i.e. a wire with in¬ 
trinsic attractive interactions (see Fig. 7a). This model 
setup can be thought of as a one-dimensional limit of the 
system analyzed in [25], where the 3D superconductor 
providing the pairing is replaced by a one-dimensional 
superconductor. 

To simplify the analysis, we assume that each wire has 
a single transverse channel. We model the system as 
two one-dimensional electron gases, labelled by z = 1,2 
(the semiconducting and superconducting wires, respec¬ 
tively). Wire 1 has repulsive short-range density-density 
interactions of strength {7 > 0, and wire 2 has attractive 
short-range interactions of strength U < 0. Both wires 
have strong Rashba-type spin-orbit interactions. As we 
will argue below, a particularly favourable case for realiz¬ 
ing the topological phase is when the spin-orbit coupling 
terms in the two wires have a similar magnitude and an 
opposite sign, such that the single-particle dispersion has 
the structure shown in Fig. 7a, and the chemical potential 
is such that the inner modes in both wires are at A: = 0. 
We will focus on this case, and comment on the effects 
of deviation from it later. We furthermore assume that 
the interactions are weak, allowing for a weak-coupling 
renormalization group analysis. 

We linearize the modes at k = 0 and at k = Ekp, de¬ 
noting the right (left) movers with spin s in wire z = 1, 2 
by R(L). = where 

^p{cr),i bosonic fields describing the 





















charge (spin) modes in wire i = 1,2. Ui^s are Klein fac¬ 
tors that impose the anti-commutation relation between 
electrons of different spins or different wires. 

In absence of tunneling between the wires, attractive 
interactions in the superconducting wire open a (non- 
topological) spin gap pinning to zero. We denote 
the magnitude of this gap by Acr, 2 - Repulsive interac¬ 
tions in the semiconducting wire give rise to a back- 
scattering term g + h.c.^ , which in terms 

of the bosonic helds can be written as 2 ^ 1 ^ cos (2(()o.,i) 
(here g > 0 is the backscattering coupling constant). As 
in the usual analysis of a repulsive one-dimensional elec¬ 
tron gas, this cosine term is marginally irrelevant and the 
spin sector in the semiconducting wire remains gapless. 

Now consider the effect of tunneling between the wires. 
Due to the spin gap in the superconducting wire, sin¬ 
gle particle tunneling between the wires is suppressed 
at energies below the spin gap in the superconductor. 
Two independent pair tunneling (Josephson) processes 
are allowed by momentum and energy conservation (see 
fig. 7b). The first process describes the tunneling of 
a pair of electrons from the fc = 0 modes of wire 1 
to the k = 0 modes of wire 2. The corresponding 
term in the Hamiltonian is given by + 

h.c.. In terms of the bosonic fields it can be writ- 
ten as ^^ 2^2 cos {26p- + (j>cr,+), where we denote 9p- = 
{0p,i - 0p,2) and cl)cr,+ = Using second- 

order perturbation theory in the single-particle tun¬ 
nelling amplitude t±, we estimate the amplitude of this 
term as gi ^ a^- 

Tunneling of a pair of electrons between the modes at 
k = ifcp of wire 1 to either the k = 0 or k = zLkp 
modes of wire 2 is given by + ^-C- or 

^21^2-t^it^U terms of the bosonic helds these 

are given by ^^ 2^2 cos {29 p- — <(>o-,±) respectively. These 
processes are suppressed with respect to the former one 
due to the mismatch in the energy dispersion for each of 
the spin havors between the wires [52], originating from 
the different spin-orbit coupling in the two wires. The 
system is therefore required to go through an interme¬ 
diate excited state at energy Agoc ~ ma^, where a is 
the spin-orbit coupling strength. The amplitude of these 
scattering processes is estimated using perturbation the- 

^2 y 

ory in t± and the interaction V as ^2 ~ —■ 

As we will see below, there is a competition between 
the two kinds of Josephson processes described above. 
The suppression of 52 with respect to gi will allow for the 
formation of a topological spin gap in the semiconducting 
wire, driving the system into the topological phase. 

We are now ready to write the low-energy effective 
action for energies below the spin gap Act ,2 in the su¬ 
perconducting wire. Neglecting fluctuations of the field 
4 >a -,2 and denoting all the quantities in the spin sector of 
the semiconducting wire simply by a, e.g. 4>cr,i —>■ 4>a, 


the bosonized Lagrangian density takes the form C = 
Ao + Aint, where 



is the quadratic part, with the index i running over the 
two wires, and 

Ant = {g cos {2(j)„) 

-I- gi cos {29p- + (j)a) + 52 cos (26'p _ - (j)a)} ■ 

(18) 


Hereafter we denote the dimensionless couplings by 
yi- We normalize the units of the imaginary time such 
that the velocity of the spin modes Ua is unity. We as¬ 
sume that the velocities of the charge modes on both 
wires are close to one and take the deviation Jrtpy = 
Up^i — 1 to be a small parameter. This allows us to rewrite 
Aq as 


T (V 0 ,) 


2 = 1,2 


Y,—Kn 


i=l,2 


27r 


id.9ppf - idr9ppy 


(19) 


where V is a 2D gradient in space-imaginary time. 

We analyze the problem using weak coupling RG. In 
the course of the RG flow additional couplings are gen¬ 
erated 

£gen = ^ (i^lV (0p,i - 0p,2) + i^2V0p,i V0p,2) 

ZTT 

( 20 ) 

We perform the RG in real space using the operator 
product expansion (OPE) formalism [47]. At each RG 
step the short distance cutoff a is increased according to 
a —>■ (1 -|- dt) a while the partition function is kept fixed 
by renormalizing the couplings. To second order in all 
couplings we obtain (see Appendix B for details): 


dK, 


-1 


dt 


y 2 I jyf + yl) 


dy , ym 

dyi . ,, K-i 

dt 

dm 

dt 




= ( 21 ) 


vm 

2 


dK, 


Pl,2 


dt 

dK, 2 

dK, 

dt 


2 I 2 

— 2/1 + 2/2 


2/2 


= -2 ivl 


■vl) 
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where ^ = {^p,i + ^p ,2 + i® the scaling dimen¬ 
sion of the scattering processes, and the velocities do not 
renormalize to this order. 

Note that yi ^2 are both relevant already to first order, 
as 2 < 2 for Ka-,Kp^i ^2 close to their non-interacting 
value of two. Assuming initially gi 3> ( 72 , in agreement 
with the discussion above, yi flows to strong coupling 
first. We denote by t* the scale at which yi becomes of 
order unity and the perturbative analysis breaks down. 
The other couplings are all marginal, but begin to flow 
significantly as t approaches t*. 

Note also the competition between yi 2 mentioned ear¬ 
lier, arising due to the term —^yy 2 {i) in the beta function 
of 2 / 1 ( 2 ) j respectively. For y > 0 (repulsive interactions in 
the quantum wire), yi tends to suppress 1/2 and change 
its sign (and vice versa). If yi and 1/2 are initially very 
different in magnitude, they may end up having an op¬ 
posite sign at scale t = t*. Their combined contribution 
enhances y while keeping its sign positive. This is exactly 
the desired situation; in order to obtain the topological 
phase y has to be positive and relevant. Having yi ^2 with 
opposite signs is equivalent to inducing a superconduct¬ 
ing gap with an opposite sign of the order parameter for 
the k = 0 modes with respect to the k = ±kF ones. For 
a fully gapped system, such a situation is exactly what 
drives the system into a (gapped) topological supercon¬ 
ducting phase in presence of time-reversal symmetry [48] . 

For t > t*, we assume that ((>1 = 26p- -|- cj)^ becomes 
strongly pinned to the minimum of the gi cosine term 
at TT. Replacing cj)i by its mean value, we obtain the 
following Lagrangian density: 

Kp,+ {V9p,+f + ^ 

{Kp,i - Alp, 2 ) (V6lp_+V()i^) +y cos {2(1)^) , 

( 22 ) 



where we ignore the contribution of the velocities Supp 
as they do not change the rest of the analysis. The Lut- 
tinger parameter of the the total charge sector is given 
by Kp,+ = j {Kp,i + Kp ,2 + 1 ^ 2 ). Note that this sector 
remains gapless as expected for a system with transla¬ 
tional invariance. The effective parameters in the spin 
sector are 


1 _ 1 
2 /= 2 /- 2 / 2 - 


1 

16 


{Kp^,+Kp^2-K2) 



(23) 


The cross term V0p^+V0o- in the Lagrangian (22) changes 
the scaling dimensions of the cosine term. However, 
treating it perturbatively we find that it does not change 
the RG equations to second order. The RG flow in the 
spin sector for t > t* then takes the standard Kosterlitz- 



FIG. 8: Phase diagram for the composite semiconducting¬ 
superconducting one-dimensional system, as function of the 
bare dimensionless pair tunneling couplings 1 / 1,2 (i = 0). Ini¬ 
tial conditions corresponding to spin-isotropic interactions in 
the quantum wire are used, | Ac (t = 0) = 1+y {t = 0), where 
the strength of the repulsive interactions in the wire is taken 
to he y {t = 0) = 0.1. We assume the non-interacting values 
for the charge sectors in each wire Kp-^ 2 (^ ~ 0) = 2. The 
scale t* at which the integration of the RG equations (21) 
is stopped and the field <j>i is replaced by its mean value was 
defined such that //i(t*) = 0.5. As argued in the text, increas¬ 
ing the spin-orbit coupling strength in the wires suppresses 1/2 
with respect to yi. As can be seen from the phase diagram 
this drives the system into the topological phase. The topo¬ 
logical phase is obtained also for j /2 yi as all the analysis 
is symmetric in 1 / 1 , 2 . 


Thouless form 


d 

dt 

dy 

dt 


{k-^) = f 

= {2 -kp 


(24) 


The topological phase is then obtained if y{t*) is posi¬ 
tive and the cosine term is relevant, i.e. if \k„ (t*) < 

^ • If fli6 cosine term is irrelevant, the spin 

sector remains gapless. 

We can now obtain the phase diagram of the system 
as a function of the bare values of the pair tunnelings 
2 / 1.2 {t = 0) as follows. The RG equations (21) are inte¬ 
grated up to the scale t* at which y^ reaches a value of 
order unity. At this scale, we replace (pi by its mean value 
and obtain the effective action (22). The flow continues 
according to Eq. 24. If the cosine term is relevant, the re¬ 
sulting phase is inferred from the sign of y{t*): y{t*) > 0 
< 0) corresponds to the topological (trivial) phase, 
respectively. 

The resulting phase diagram is shown in Fig. 8. We 
used initial conditions that correspond to spin-isotropic 
interactions in the quantum wire, \Krj (t = 0) = 1-1- 
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y {t = 0). We set y [t = 0) to 0.1. For simplicity we as¬ 
sume the non-interacting values for the charge sectors in 
each wire (^ = 0) = 2. The scale t* was defined 
such that yi{t*) = 0.5. Neither of these choices changes 
the phase diagram qualitatively. 

We find that the cosine term always flows to strong 
coupling opening a spin gap in the semiconducting wire. 
The topological phase is obtained for yi sufhciently larger 
than corresponding to pair tunneling of the outer 
modes g 2 being suppressed with respect to the pair tun¬ 
neling of the inner modes gi. In our microscopic system 
this situation can be achieved due to the form of the spin- 
orbit coupling in the wires. Setting the chemical poten¬ 
tial to the point indicated in Fig. 7b is optimal for this 
purpose; large enough spin-orbit coupling will eventually 
drive the system into the topological phase, since it sup¬ 
presses 52 but not gi. Note that the topological phase 
can also be obtained for 1/2 yi, as all the analysis is 
symmetric in yi, 2 - 

V. DISCUSSION 

We have shown that a gapless topological phase, 
with exponentially localized edge states, can exist in a 
strictly one-dimensional time-reversal symmetric (class 
Dill) system with a conserved number of particles. The 
localization of the edge states arises due to a gap to sin¬ 
gle particle excitations in the bulk of the system. The 
edge states are characterized by an anomalous relation 
between the fermion parity and time reversal symmetry: 
time reversal flips the local fermion parity at the two 
edges, just as in the case of a Dill mean-field topological 
superconductor [16]. 

It is interesting to note that gapless topologically pro¬ 
tected phases are possible in other symmetry classes in 
one dimension, as well (see, e.g.. Refs. [49, 50]). In par¬ 
ticular, class BDI (with particle-hole symmetry and time 
reversal that squares to -1-1) is likely to support similar 
phases, although there are less distinct phases than in 
the non-interacting case [51]. As we have shown here, as 
long as there is a gap for single fermion excitations in the 
bulk even in the absence of long-range order, the corre¬ 
sponding topological phase may be robust in a gapless 
one-dimensional system. 

We discussed a possible realization of the phase by 
proximity coupling a semiconducting wire with spin-orbit 
coupling to a superconducting wire. Another realization 
of the phase may be possible in quasi-one-dimensional or¬ 
ganic conductors where spin anisotropic interactions are 
believed to be present. As was analyzed in Ref. [34], these 
interactions could also drive the system into a phase with 
dominant triplet superconducting correlations, that we 
identified as a topological phase. In such a system, inter¬ 
chain hopping eventually drives the system into a three- 
dimensional long-range ordered superconducting phase; 


however, if the system is very anisotropic, the properties 
of the topological phase (e.g., the presence of zero-energy 
surface states) may already be apparent at intermediate 
temperatures above the critical temperature, where the 
system is described as nearly-decoupled fluctuating one¬ 
dimensional superconductors. 

It is important to note that, unlike the fully gapped 
topological superconducting phase which is immune to 
non-magnetic disorder (as long as its strength does not 
exceed the superconducting gap), the gapless phase can 
be sensitive even to weak disorder. In the topologi¬ 
cal gapless phase 2kp density fluctuations are gapped. 
Therefore the dominant effect of an impurity potential is 
due to the Akp density fluctuations term. The effect of 
this term depends strongly on the nature of interactions 
in the system, parametrized by the Luttinger parame¬ 
ter of the charge sector, Kp. For Kp sufficiently large 
(in the conventions used here, Kp > 1, where Kp = 2 
corresponds to a non-interacting system), backscattering 
due to a weak point-like impurity potential flows to zero 
under RG, leaving the topological phase intact. If the 
interactions are too strongly repulsive, corresponding to 
Kp < 1, any weak impurity potential becomes relevant. 
In the presence of many such impurities, the system be¬ 
comes localized; this destroys the topological nature of 
the phase, and gaps out the edge states. 

We conclude that the topological phase is robust in 
the presence of weak disorder, as long as the repulsive 
interactions in the system are not too strong. A metallic 
gate placed near the system may be used to screen the 
long-range part of the Coulomb interactions, thus making 
the topological phase more stable to disorder. 
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Appendix A: Tunneling density of states at an edge 
of a gapless topological system 

In this section we describe the calculation of the matrix 
element given in Eq. 9, which contributes to the a; —>■ 0 
tunneling density of states (TDOS) for a system in the 
gapless topological phase. 

Consider a topological region that extends from x = 0 
to X = L, with a trivial region on its either side. The 


11 


bosonized Hamiltonian describing the system is given by 

+ J dx^^^^cos{2(l)a), (Al) 

where g (a:) = 5 > 0 in the topological region 0 < x < L, 
and g (x) = go < 0 in the trivial regions, x < 0 and 
X > L. We assume the couplings g,go are large enough 
such that (f>a- is pinned close to the minimum of the poten¬ 
tial in each region. This allows us to expand the cosine 
around nni^2 in the trivial regions and around tt (to + 5 ) 
in the topological region, where ni 2 ,m are integers. We 
consider a variational Hamiltonian where the cosine term 
is replaced with a quadratic potential, consistent with 
such an expansion. There are two distinct configurations 
of the potential, corresponding to different states: 


where an expectation value with respect to the ground 
state of Hq is assumed in the final expression. Writing 
the single particle creation operator as 

4-1 (x) ~ 

r 

where (x) — Op {x) + \r<p^ {x) — 9^ (x) and 

the sum over r = ±1 stands for right and left movers 
respectively, the matrix element becomes 

r 

To proceed we therefore need to diagonalize and find 
the expansion of the bosonic fields pp^a and Op^a in terms 
of its eigenmodes. Note that since the spin and charge 
sectors are decoupled, the expectation value above can 
be written as a product 



Vo{x) 


1 {-2go(j)l X < 0, x> L 

27r2a2 (25 - f 0 < a: < L 




n 

a—p,<7 




and 

^ f-2go (0a - tt)^ a: < 0 
~ 27^2^2 I ‘^0 (0a — 2 ) 0 < X < L ■ 

[-2go0a X>L 

The respective variational Hamiltonians are denoted by 
Hq^i and their ground states by |0) and |1). Saddle point 
configurations of pa in these states, 0a, 0,1 are de¬ 

picted in Fig. 1. We denote the difference between them 
by A0cr = pa,i — 0a, 0 - Note that the energy of the spin 
sector for these two configurations is the same, and they 
are therefore nearly degenerate. The difference in energy 
arrises due to the charging energy (which scales as 
tending to zero for large enough system), as in the state 
|1) the total spin in the system is corresponding to an 
odd number of particles, while in the state |0) the total 
number of particles is even. Since the density of the par¬ 
ticles in the system is given by —^dxpp, the respective 
saddle point configurations of 0p in these two states dif¬ 
fer by a constant gradient term A0p = tt (l — j) in the 
region 0 < x < L. The two ground states are therefore 
related by a unitary transformation |1) = |0), where 

6 = da and 

0!. — p,(7 

1 

Sa = — dx' Apa (x) dx'Oa {x) ■ 

^ j — 00 

To see this note that under this transformation the field 

pa ^ pa + Apa- 

The matrix element that is expected to give the largest 
contribution to the TDOS at w —>■ 0 can now be written 
as 

(cr)|0) = (e%| (x)), 


where {■■■)a=piy denotes the expectation value in the 
ground state of the charge and spin sectors respectively. 

We begin by diagonalizing the spin sector, i?o,a- For 
simplicity we take the spin gap in the trivial regions to 
be inhnite, |go| 00 (equivalently we can assume that 
the topological region is surrounded by vacuum). This 
pins the field pa and the current dxOa at the boundary 
to zero, i.e. Pa\x=Q l ~ ® d)xOa\x=Q l = 0. A mode 
expansion for the helds pa , da then takes the following 
form: 



where the expansion for pa is around the constant saddle 
point solution pafi{x), and are bosonic creation 

and annihilation operators satisfying the commutation 


relations 


= 0 


k.k' • 


The Hamiltonian translates into 


Ho,a = 


I 00 

X ( 


^4 + 






(4 + 4 ^) 


where Ak = Ua^ + 


\g\K,, 




and a Bo- 


Tk ’ ^ Tra^ tt/c ’ 

goliubov transformation = oikb^ + I3kb\ can be used 
to bring it into a diagonal form 


Ho,a = ^kdlK + const. 


The coefficients in the transformation are ak = 
^ 5 ( 1 ^ + ^)’ l)- The eigen-energies 
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are given by Ek = y/— Bf and 6^, 6^, are the bosonic 
creation and annihilation operators of the corresponding 
eigen-modes. 

The helds (j)a, are given by linear combinations of 
these creation and annihilation operators. Hence, the 
expectation value can be calculated using the identity 


Denote the spin gap in the system by Ao- 


2\g\u„K„ 


and the inverse correlation length by k = 

In the limit j^ol —>■ oo the difference between the saddle 
point configurations is 


The absolute value of the matrix element squared is 
therefore 



where we denote a = —'r \Kp and ^ is the 

correlation length. 


Appendix B: RG flow equations obtained using the 
operator product expansion 


{ TT X < 0 

’r -'S(hr” 

0 X > L 

We can then write 6cr = Qr, (0) -I- where A0o- = 
dxe~‘^^dx0a (x) for kL ::$> 1. Note also that we are in¬ 
terested only in the absolute value of the matrix element. 
The expectation values {(pa-(x) 9a (x)) and (^pa (x) Sa'j 
are imaginary and only contribute a constant phase shift 
to the oscillations in 2kF of the matrix element squared. 
The asymptotic behavior is dictated by the function 


Fa=(^ (9a (x) -Sa'j ^ -h ^ {pi {x)) . 



In this appendix we give further details on the deriva¬ 
tion of the RG flow equations (21). As mentioned in the 
main text, we perform the RG in real space. At each 
step we rescale the short distance cutoff a, defined as the 
minimal distance between two operators in the theory, 
according to a —?> (1 -I- dt) a. We use the operator prod¬ 
uct expansion (OPE) to replace pairs of operators which 
are within the new short distance cutoff, allowing us to 
recover the original action with renormalized couplings. 

The general form for an OPE of two operators Oi, Oj 
is 

: O, (ri) :: Oj (ra) ^Cijk (n - ra) : Ok b 

k b / 

where : Oi : denotes normal ordering of the operator. 
The equality is valid only when both sides of the equa¬ 
tion are considered inside a correlation function with an¬ 
other operator (or set of operators) at a distance much 
greater than \ri — ra| from the operators Oij . The func¬ 
tions Cijk (ri - ra) have the form Ci^k = ^ 

where dij^k are the scaling dimensions of the correspond¬ 
ing operators and Cijk are the OPE coefficients which are 
pure numbers. The beta functions to second order of the 
respective couplings gij^k are then given by [47] 


Taking the continuum limit 0, the sum over k can 
be written as an integral over q = ^. Performing the 
integral and considering the asymptotics for ^ a: A 
we find that 


Fa 


TT 1 


KX. 


A similar calculation for the charge sector gives the 
following asymptotic behavior ior a x L 


dgk 

dt 


{d dk) gk ^ ^^ Ojkgigj: 
i,3 


where d is the dimension of the problem. 

We demonstrate the OPE explicitly for the pair of op¬ 
erators Oi = \79p-\7pa and O 2 = cos {29p- pa), both 
of which we treat perturbatively. We start by writing the 
cosine as a sum of exponents, and expending the expo¬ 
nent as a series 




(9p{x)-5p^ \ + i (<(.2 (x)) 


X + ^Kg 1 In ( - ) -p ^ In 


2 \Kg 


(D 


Kg 
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The normal ordered product of the operators : Oi (ri) :: and 

O 2 (^ 2 ) : will then contain a summation over terms of the 
form 


: (n) (n) :: {r^) C”™ (^2) : • 

In each such term, we can contract diOp^- (ri) with one of 
the (r 2 ), giving m possible configurations. For each 
of them, there are n — m ways to contract dicj)^ (ri) with 
one of the (pa- (^ 2 ). These contractions are given by 

{d.PAri)ct>a (r2))=-JiF.f^^^ 

^ \ri-r2\ 


{di0p- (ri)6»p_ (r2)) = 

{di 9 p^i (ri) 6»p_i (r2)) + (9i0p,2 (ri) 6»p,2 {r2)} = 

« l^pl + ^p2 ) I ,2 

^ |ri-r2| 

respectively. The product of the operators now becomes 


VOp^.VPa := ^Ka {K;} + K;]) 


I 


ki - r2k^^!' 




Rewriting the sum over m as 

: (20p,_r-:= 

m ^ '' 

2n (n-l) : (26»p _ + (j)aT~'^ ■ 

and the sum over n as 

,_. An —2 

n ^ '' 

we obtain 

: VOp-Vcpa ■■ cos (20p,_ + pa) := 

— {^p,i + '^p,2) "i • *^0® (2^p,- + ko-) b 

^ ki-r2| 

where we identify —^Ka {K~l + K~l) as the OPE coef¬ 
ficient. This results in the following contribution to the 
beta function of j/i 

^ = \Ka{K;}+K-l)Km- 

To second order in the coupling constants, we can use 
the non-interacting values for the Luttinger parameters 
Ka = Kp^i ^2 = 2 and write 
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